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CHAPTER  I 


INTRODUCTION 

The  shallow  arch  represents  one  of  the  simplest  realistic  structures 
from  which  many  of  the  features  of  elastic  instability  theory  can  be 
illustrated.  As  a  result  of  this  simplicity  the  various  existing  theories 
which  have  been  proposed  to  describe  the  behavior  of  an  arch,  as  well  as 
the  methods  of  constructing  such  theories,  have  received  considerable 
attention  for  over  fifty  years.  More  recently,  due  to  the  development  of 
new  types  of  materials  being  considered  for  use  in  submersible  and  space- 
type  structures,  it  has  become  desirable  to  include  anisotropic  material 
properties  in  the  theory. 

The  current  investigation  will  be  devoted  primarily  to  the  problem 
of  structural  behavior  of  an  elastic  orthotropic  shallow  arch  including 
transverse  shear  and  transverse  normal  stress  effects  and  undergoing 
finite  deflections.  The  effect  of  initial  imperfections  is  also  considered 
in  order  to  investigate  possible  bifurcation  of  the  arch. 

The  general  theory  derived  here,  which  considers  transverse  shear  and 
transverse  normal  stress  existing  in  an  orthotropic  arch  is  later  specialized 
to  the  case  of  an  isotropic  arch  with  a  consideration  of  shear  deformation 
and  is  further  specialized  to  the  still  simpler  theory  of  the  classical  arch 
which  is  then  compared  to  a  theory  derived  independently.  For  the  case  of 
asymmetric  deformation  the  transition  from  the  initial  state  to  the  buckled 
state  is  shown  to  be  connected  with  geometrical  imperfections  present  in  the 
structures. 

A  survey  of  the  literature  on  arches  reveals  that  many  papers  have  been 
published  pertinent  to  the  classical  theory  of  arches.  A  general  discussion 
and  significant  references  up  to  1955  covering  buckling  problems  of  circular 
rings  and  arches  is  given  by  S.  Timoshenko  and  J.  M.  Gere  [1],  These  authors 
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shew  that  in  the  case  of  very  flat  arches,  buckling  in  which  axial  strain  is 
considered  may  occur  at  a  smaller  load  than  inextensional  buckling.  Timoshenko 
assumed  that  the  center  line  of  the  deflected  arch  as  well  as  the  initial  no- 
load  shape  is  a  half  wave  of  the  sine  curve  and  arrived  at  a  very  simple  solu¬ 
tion.  A  very  complete  investigation  of  structural  behavior  of  a  shallow  arch 
is  attributed  to  Y.  C.  Fung  and  A.  Kaplan  [2].  These  authors  presented  a 
series  solution  for  the  case  of  an  arch  having  a  sinusoidal  configuration  at 
no-load.  They  also  found  that  instability  m3y  be  associated  with  a  bifurcation 
into  an  asymmetric  buckling  mode.  The  same  approach  was  employed  earlier  by 
C.  Biezeno  [3]  but  for  the  case  of  a  circular  arch  loaded  by  a  concentrated 
normal  force  at  the  center.  K.  Federhofer  solved  the  buckling  problem  of  a 
parabolic  arch  in  1934  [4]  and  discussed  fhe  dynamic  problems  of  arches  and 
rings  later  [5].  K.  0.  Friedrichs  [6]  and  K.  Marguerre  [7]  derived  the 
equations  for  a  circular  arch  under  uniformly  distributed  pressure  by  varia¬ 
tional  principles.  K.  0.  Friedrichs  pointed  out  that  the  uniform  load  corre¬ 
sponding  to  unsymmetric  buckling  could  be  much  lower  than  that  for  symmetric 
buckling  «nd  in  addition  advanced  a  criterion  for  the  buckling  load  as  "the 
lowest  load  at  which  a  buckled  state  exists  having  the  same  energy  level  as 
the  unbuckled  state  at  the  same  load."  The  so-called  equal-energy  criterion 
due  to  A.  Gjelsvick  and  S.  R.  Bodner  [8],  although  net  generally  accepted  as 


a  proper  buckling  criterion  may  lead  to  an  "energy  load"  that  represents  a 
significant  lower  bound  to  the  buckling  load.  h.  Schreyer  and  E.  Masur  [9] 
gave  the  solution  for  the  structural  behavior  of  a  clamped  circular,  shallow 
arch  together  with  a  detailed  analysis  of  the  various  buckling  criteria  by 
considering  the  buckling  equations  and  using  the  change  in  potential  energy 
between  the  buckled  and  unbuckled  states. 

Notable  contributions  to  the  theories  of  higher  oroers  of  accuracy  which 
include  the  effects  of  either  or  both  transverse  normal  stress  and  shear 
deformation  are  due  to  F.  Hildebrand,  E.  Reissner  and  6.  Thomas  [10],  A.  E. 
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Green  and  W.  Zerna  [11],  E.  Reissnar  [12],  and  by  P.  M.  Naghdl  [13].  In 
reference  [10],  orthotropic  shells  are  considered  and  the  various  theories 
which  are  derived  from  different  approximations  in  the  assumptions  are  com¬ 
pared  and  discussed.  A.  E.  Green  and  W.  Zerna  concluded  that  the  inclusion 
of  shearing  forces  is  likely  to  be  of  importance  mainly  in  the  case  of 
edge  effects.  E.  Reissner's  work  on  obtaining  stress-strain  relations  by 
removing  the  first  approximation  of  Love  [14]  results  in  some  simplifications 
as  compared  with  those  given  in  [10]  and  [11].  Later  P.  M.  flaqhdi  included, 
in  addition  to  the  transverse  shear  and  transverse  normal  stress,  the  rota¬ 
tory  inertia  in  isotropic  shell  analysis.  E.  Reissner  [15]  in  a  separate 
work  discussed  a  problem  concerning  the  torsion  of  a  rectangular  plate  and 
this  work  has  become  the  standard  of  comparison  in  the  linear  bending  theory 
of  elastic  flat  plates  which  includes  transverse  shear  effects.  In  the  large 
deformation  range,  P.  Wilson  [16]  investigated  the  elastic  circular  flat  plate 
with  the  transverse  shear  effect  taken  into  account  and  R.  R.  Archer  [17] 
derived  a  nonlinear  shear  deformation  theory  for  isotropic  shells.  In  the 
present  investigation,  the  initial  imperfections  are  also  introduced  to 
achieve  better  correspondence  with  the  true  structure.  This  leads  to  lower 
buckling  loads. 
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CHAPTER  II 
CLASSICAL  THEORY 

I 1.1.  Introduction 

In  this  chapter  the  classical  theory  will  be  developed  and  the  critical 
loads  corresponding  to  various  values  of  the  arch  parameters,  X  ,  will  be 

I 

determined.  The  von  Karman  type  nonlinear  strain-displacement  relations  are 
used.  The  arch  considered  here  is  circular,  thin  and  shallow.  In  order  to 
obtain  a  theory  for  the  arch  which  is  similar  to  the  Euler-Bernoulli  beam 
theorem,  we  make  the  following  assumptions:* 

(a)  A  normal  to  the  undeformed  arch  middle  surface  is  deformed  without 
change  of  length  into  a  normal  to  the  deformed  middle  surface. 

The  change  in  direction  of  the  normal  to  the  middle  surface  is 
neglected.  Thus  the  change  in  direction  of  the  transverse  load 
is  then  neglected. 

(b)  The  stress  in  direction  of  the  normal  to  the  middle  surface  is 
small  compared  to  the  other  stresses. 

(c)  The  ratio  of  thickness  h  to  the  radius  of  curvature  R,  is  negligibly 
small  compared  to  unity. 

(d)  The  magnitude  of  transverse  deflection  is  of  the  same  order  as  the 
thickness  of  the  arch. 

(e)  The  load  is  considered  to  be  applied  at  the  middle  surface. 

The  theory  is  called  classical  in  the  sense  that  it  is  based  on  these  assump¬ 
tions  which  result  in  ignoring  the  effects  of  shear  and  the  transverse 
normal  stress. 

Based  on  these  assumptions,  the  governing  differential  equations  are 
derived  by  the  principle  of  minimum  potential  energy.  A  finite  difference 

*  This  assumption  is  similar  to  the  Kirchhoff  assumption  of  plate  theory  or 
Love's  first  approximation  of  shell  theory. 
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scheme  together  with  Newton's  iteration  technique  is  used  to  solve  the  equa¬ 
tions.  Initial  imperfections  are  also  considered. 


II. 2.  Strain-Displacement,  Stress-Strain  Relations 

Let  the  origin  of  an  x-z  coordinate  system  be  located  at  one  end  of  the 
arch  with  x  being  the  coordinate  extending  along  the  middle  surface  cf  the 
arch  and  z  being  the  distance  of  an  arbitrary  point  measured  along  the  outward 
normal  to  the  middle  axis  of  the  arch  (see  Fig.  1(a)).  For  strain-displace¬ 
ment  relations,  let  us  employ  the  von  Karman  type  approximation  in  which  the 
only  non-linear  terms  retained  are  those  corresponding  to  quadratic  effects 
in  the  transverse  displacement  and  its  derivatives.  Thus  the  strain-displace¬ 
ment  relation  is 


e 


x 


1  r  3W  -.2 

2r  3X  J 


1/3WOn2 


2'  3  X 


-r 


(2.1) 


where  U,W  are  the  displacement  components  at  an  arbitrary  point  in  the  arch  in 
the  x,z  directions  respectively  and  w  is  the  initial  imperfection  in  the  radial 
direction  (see  [18]). 

The  Euler-Bernoulli  type  theorem  as  stated  in  section  2.1. leads  to 


U(x,z)  =  u(x)  -  z[w( x)  +  wQ(x)],x 
H(x)  =  w(x)  +  wQ(x) 


(2.2) 


where  u(x)  and  w(x)  are  the  components  of  displacements  at  the  middle  surface. 

By  substituting  (2.2)  into  (2.1),  the  strain-displacement  relation 
becomes 


Ev  =  u,v  +  j£  +  \  w,?  +  wft  „  w,v  -  z  (w,vv  +  w„  VJ 


X  *x 
■  8x  +  zKx 


X  o,x  X 


'XX  0,XX' 


(2.3) 


where  e 


wo.xw'x  and  Kx  *  -w*xx  -  Vxx  are  018  m1ddle 


surface  strain  and  the  change  of  curvature,  respectively. 


(a)  COORDINATE  SYSTEM  a  GEOMETRY 


(b)  ARCH  ELEMENT 


FIG.  I  THE  CONFIGURATION  OF  THE  ARCH 
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The  stress-strain  relationship,  assuming  Hooke's  law  holds,  is 


a  =  Ee 
X  x 


(2.4) 


or  conversely 


ex  =  E  * 


(2.5) 


I I. 3.  Derivation  of  Equilibrium  Equations  and  Natural  Boundary  Condi lions 
The  strain  energy  density  UQ  is  a  function  of  the  strain  component  [19], 


i  .e.: 


U0  ■  U0(ex)  (2.6) 

If  the  displacement  vector  (u,w)  receives  variations  (6u,  6w),  the  strain 
tensor  assumes  a  variation  6ex  and  the  corresponding  variation  of  UQ  is 

3U 

iUc  ’  *£  5V 

It  has  been  shown  that  the  following  relation  holds  for  a  Hookean  material  in 
the  large  deformation  range  [19], 


°x  3e. 


Therefore 


SU  =  a  <5e 
o  xx 


(2.7) 


Substituting  (2.3)  into  (2.7),  one  finds 


sUo  ■  °x:,'Ex  +  ZSKx!- 


Let  the  first  variation  of  the  strain  energy  of  the  entire  arch  be  6Uy.  Then 

L  h 

aij_  =  f  f  F  ah  H? dx 
T  Jn  J ± 


6UQdzdx 


,  h 

rL  r  7? 


=  |  J  «*x(«c°  +  z6Kx)dzdx  (2.8) 
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or 


L  h 

{UT  =  J J  h  e[4  <  +  z(cx  5Kx  +  Kx  ^  *  z2rx  sKx]dzdx 
0  ~2 

Upon  integration  with  respect  to  z,  one  finds 
SUT=  £  (Eh  'x5ex  +  TTkx  5Kx>dz 


(2.9) 


■  iUn,  +  sUb 


L. 

where  6Um  =  j  Eh  e°  6e°  dx  is  the  membrane  strain  energy  and 

L  t 
f  Fn 

b  =  J  TF  6Kx  bending  strain  energy. 


The  force  and  moment  resultants,  fl  and  M  (per  unit  width  of  arch) 
shown  in  Fig.  1(b)  take  the  form 


h 


N 


M 


=  u  c..dz 

)  II  A 

2 

(2.10) 

=  f  h  !t,d!' 

(2.11) 

"2 

which  upon  integration  with  respect  to  z  become 


n  =  Eh(u’x  +  l+  wo,x  w’x  +  Iw’x}  =  Eh  ex 

o  3  (2.12) 

m  -  4.  ^  it 

M  -  +  w),xx  -  72- Kx 

Substituting  these  force  and  moment  resultants  into  the  first  variation  of 
the  total  strain  energy  (2.9),  we  get 

6UT  =  L  CM  6ex  +  H-y  ^ 


(2.13) 


Now,  substituting  (2.3)  into  (2.13),  we  get 
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-fO(w,xx)]}  dx 


(2.14) 


The  potential  energy  due  to  the  external  load  q  is 

n  =  -  j  q(w  +  wQ)  dx 

and  the  total  potential  energy  n  is  then 
n  =  UT  +  n. 


By  variation  of  n,  we  get  the  equilibrium  equations  and  the  associated 
natural  boundary  conditions.  The  equilibrium  equations  are: 

N,x  =  0  (2.16) 

“  w’xx  "  wo,xx^  "  H’xx  “  q  ~  0 


and  the  natural  boundary  conditions  are: 
fl  =  0  or  u  prescribed  at  x  =  0  and  L 
•lw,  +  M,  =  0  or  w  prescribed  at  x  =  0  and  L* 

A  } \ 

M  =  0  or  w,  prescribed  at  x  =  0  and  L 

A 


r .  .i  .xi  -  /  n  1  rt\  j..  x.  fry  i£\  ... 

oUuboi  outing  \c.\c)  iiilu  We  Qeo 


U,XX  +  R  w,x  +  wo,x  w,xx  +  wo,xx  w’x  +  w*x  w’xx  ”  0 


1 2 ^wo  +  w^’xxxx  *  (U’x  +  R  +  wo,x  w’x  +  2  w,x^ 


1  ...  2, 


(S  -  wo,xx  -  w'xx>  -  1  *  0 


(2.17) 


II. 4.  Solution  to  the  Equations 

Let  us  introduce  the  following  nondimensional  variables 


*  Nw,  +  M,  corresponds  to  shear  forces  V  in  Fig.  1(b). 

A  A 


~~9 - 


/V 


and  X  = 


The  equations  (2.17)  then  take  the  form: 


Ss  +  A“’s  +  wo,s  “’ss  +  wo,ss  °’s  +  °'s  “'ss  ■  0 


(2.18a) 


iVa  +  w0)*ssss  +  (x  “  ”o,ss  ■  a‘ss)  ^»s  *  Xa  *  v/o,s  a’s  *  “  Q  =  °*  (2*18b) 

where  Q  =  qL4/Eh4 

Let  us  next  linearize  these  equations  in  order  to  obtain  a  system  that 
can  be  treated  by  Newton's  iteration  scheme  [20].  Regarding  the  displacement 
components  $  and  a  and  their  derivatives  as  independent  variables,  let  us 
expand  the  nonlinear  terms  of  those  variables  Into  Taylor  series  using  as  a 
starting  point  the  no-load  equilibrium  configuration  and  retaining  only  the 
nonlinear  parts.  Typical  expansions  of  the  terms  are  as  follows: 


°’s  a’ss  =  ®»s  “’ss  +  “’sKs  “  «.*,)  +  “’ss(ct»s  “  “>s> 


=  a,,.  a,cc  +  o,cc  a,c  -  a>e  a,pt. 


a’ss  a’s  =  2“’ss  --S  a’s  *  ~’s  a’ss  ‘  2“*ss  ®»s 


where  the  double  barred  quantities  are  the  estimates  of  the  same  unbarred 
quantities.  The  linearized  versions  of  the  equations  (2.18)  are 

+  Xa,_  +  W„  _  a,__  +  W.  a.  +  a.  _  a,„  +  a,,  a._  =  a.  a.  „ 

T  ss  ’s  o,s  ss  o,ss  s  ss  s  ’s  SS  S  S£ 

(2.19a) 


12  ’ssss 


+  (X  -  a.  c  -  Wn  __)  ($,-  +  Xa  v  W„  c  a,  +  a.  a,J 


’SS  0,SS'  'T’S 


0,S  S  ’S  rSJ 


-  1  ? 

*  In  [6],  because  of  (2.16),  (q>,s  +  Xo  wQ,s  o,s  +  j  a,p  was  treated  as  constant. 


n 


■<“-s  +  J"  +  \  "4 )  “-ss  =  '  T7  Ko,ssss  +  I  x“4  -  *,  "-ss 


x”-ss  ■  ”>ss  "4  •  7wo,ss  "4  •  wo,s  ”’s  ~'ss  +  Q 


(2.19b) 


The  quantities  without  bars  in  (2.19)  are  approximated  by  central  differences. 
If  we  select  N  internal  grid  points  in  the  x  direction,  the  grid  increment  is 
then  1/N+l  and  equation  (2.19)  becomes 


22  n  2—2  n  —  n ■■ 

D  6  6.  +  XiS  &•  +  D  W  S  cu  +  t?  W  5  ctj  +  7r  Ot,  6  a* 

S  V1  2  S  I  0,S  S  1  2  O.SS  S  ,  2  JSS  S  1 


,  n2  -  .2  _  -  - 

+  U  a,  o  a •  =  cu  cu  . 
*S  S  1  S  ‘ss 


(2.20a) 


D4  .4 


I'sV1*-  “o.ss  •  “*ss>  (f  ss  *1  +  x“i  +  I  «o,s  5s  °i  +  Fs  {s 

-«?•,  +  x“  +  ?  “4>  1)2  4  “1  =  '  17  “o.ssss  +  7  x"4  '  N  "'ss 
-  x"  «.ss  -  “,ss  -  Rss  “4  -  50,S  “’s  "'ss  +  Q  (2-20b) 


where  D  =  N  +  1 


6s  *i  =  ^1+1  “  ^i-1 

°'s  *1  =  ^i+l  "  2*i  +  ^i-l 

6s  *1  =  *1+2  "  4*i+l  +  6*i  "  4*i-l  +  *i-2* 


Equation  (2.20)  represents  a  system  of  2  N  simultanaous  nonlinear  algebraic 
equations  which  are  solved  by  Newton's  iteration  scheme  osing  a  digital 
computer.  Ordinarily  this  has  been  done  by  input  load  increments  as  the 
known  quantity  to  solve  for  $  and  c.  However  the  procedure  fails  to  give 
any  information  on  the  load-deflection  relations  after  the  first  peak  load 
is  reached,  since  the  slope  is  zero  at  this  point  and  the  computer  yields 
unstable  solutions.  For  the  present  case  involving  a  consideration  of  imper¬ 
fections,  the  extension  of  the  load-deflection  curve  past  its  peak  becomes 
desirable  to  be  able  to  obtain  the  lower  buckling  load. 
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If  deflection  instead  of  load  is  employed  as  the  input  and  load  regarded 
as  the  unknown  quantity,  the  above-mentioned  difficulty  is  avoided  and  the 
load-deflection  curve  can  be  extended.  But  since  the  buckled  mode  shape  of 
the  arch  is  yet  to  be  determined,  the  input  of  deflections  is  not  realistic. 

As  an  alternate  approach  the  summation  of  all  radial  displacements,  i.e.,  the 
change  of  area  between  the  arch  and  its  base  line  can  be  employed  by  introduc¬ 
ing  an  additional  equation  relating  the  area  to  the  radial  deflections,  namely: 


AREA 


1  ? 
—  I  a, 
N+l  i=l 


This  equation  together  with  the  previous  2N  equations  forms  a  system  which 
has  the  right  number  of  equations  after  the  idditional  unknown  Q  is  added. 
Let  us  write  the  system  in  matrix  form 


[A]  [X]  =  [B] 

where  [A]  is  a  coefficient  matrix,  [B]  is  a  constant  column  matrix  and  [X] 
the  unknown  column  matrix.  The  loading  term  Q  in  (2.20b)is  moved  to  the  left 
hand  side  of  the  equation  and  is  included  in  the  unknown  matrix  [X].  The 
double  barred  quantities  in  [A]  and  [B]  are  assumed  to  be  zero  at  the  beginning 
of  the  iteration  or  more  precisely,  the  displacement  components  are  assumed  to 
be  zero  at  the  beginning  of  the  iteration.  Then  the  independent  area  para¬ 
meter  AREA  in  [B]  is  given  successive  increments  and  the  displacements 
determined.  If  the  :alculation  for  displacements  converges,  a  further  incre¬ 
ment  is  employed. 

If  we  neg'pet  the  Imperfection  terms  in  the  equations,  the  solution 
for  radial  displacement  u  in  general  is  symmetric  with  respect  to  the  mid¬ 


point  of  the  arch,  since  the  equations  are  expressed  in  central  differences 
and  the  boundary  conditions  are  symmetric.  However  the  symmetric  characteris¬ 
tics  of  the  syscem  are  of  course  disturbed  if  we  impose  an  antisymmetric  mode 


of  initial  no-load  imperfection.  The  resulting  solution  is  then  an  asym¬ 
metric  one  and  instability  occurs  at  a  lower  load  than  in  the  case  of 
symmetric  buck! ir. j  of  an  initially  perfect  arch. 

For  the  hinged  end  arch,  an  antisymmetric  sine  curve  is  assumed  to 
represent  the  no-load  configuration,  viz: 

*o  =  *oo  s1n  2l,s 

where  wQ0  is  the  amplitude  and  a  numerical  value  of  0.04  is  used  for  numeri¬ 
cal  analysis.  For  the  clamped  end  arch  a  polynomial  corresponding  to 
clamped  end  conditions  and  antisymmetric  with  respect  to  the  mid-point  of 
arch  is  used.  The  polynomial  is  wQ  =  wQ0  [As  +  Bs  +  Cs  +  Ds  +  Es  ], 
where  A  =  2.048910  x  10~8,  B  =  -  1.137778  x  102,  C  =  2.844444  x  102, 

D  =  -  2.275556  x  102,  E  =  5.688889  x  10  and  w  „  =  0.04. 

oo 


I I. 5.  Discussion  of  Results 


Typical  load-deflection  curves  for  hanged  end  and  clamped  end  arches  are 
shown  in  Figs.  2  and  3  respectively.  There,  the  asymmertic  solution  is 
obtained  b j  imposing  antisymmetric  initial  imperfections  while  the  symmetric 
solution  is  the  one  obtained  by  setting  the  amplitude  w0Q  equal  to  zero.  To 
determine  the  load-deflection  relations  shortly  before  point  "C"  of  Fig.  2 


has  been  reached,  the  imperfection  wQQ  (introduced  to  force  the  asymmetric 
form  to  appear)  is  taken  to  be  zero  and  the  return  branch  cb  is  found.  The 


bifurcation  point  "b"  where  the  returning  curve  meets  the  symmetric  curve  is 


about  25  percent  lower  than  the  point  "a"  in  this  case. 


The  variation  of  load-deflection  relation  with  the  geometric  parameter 
X  is  shown  in  Figs.  4  and  5.  Figs.  6  and  7  show  how  the  critical  load 
increases  with  a.  From  these  figures  the  following  conclusions*  can  be  drawn: 
(1)  For  the  hinged  end  arch: 

(a;  For  A  s  6,  no  buckling  occurs; 


*  figures  are  approximate  values. 


IRtSJ  Kflasu  ramo 


FIG.  4  LOAD-DEFLECTION  CURVES  FOR  THE  HINGED  ENDS 


SL 


tc)  x=!8 

FIG.  5  LOAD-DEFLECT!ON  CURVES  FOR  THE  CLAMPED  END 


CURVES  FOR  SYM.  AND  ASYM.  SOLUTIONS 
CLAMPED  ENDS 
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(b)  For  6  <  A  s  10.8  (10.5  by  Y.  C.  Fung  and  A.  Kaplan  [2]), 
the  buckling  mode  is  symmetrical; 

(c)  For  A  >  10.8,  the  asymmetric  buckling  criterion  holds.  The 
asymmetric  solution  found  here  agrees  very  closely  with  Y.  C. 
Fung  and  A.  Kaplan's  solution  [23  which  is  also  plotted  on 
Fig,  6. 

(2)  For  the  clamped  end  arch: 

(a)  For  A  <  10.0  (11.4  by  H.  Schreyer  and  E.  Masur  [9]),  no 
buckling  occurs; 

(b)  For  10,0  s  A  <  25  (22.96  by  H.  Schreyer  and  E.  Masur  [9]), 
the  arch  behavior  is  governed  by  the  synmetric  criterion; 

(c)  For  A  2  25,  the  asynraetric  buckling  criterion  holds. 
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CHAPTER  III 

TRANSVERSE  SHEAR  DEFORMATION  AND  TRANSVERSE 
NORMAL  STRESS  THEORY  TOGETHER  WITH  A 
CONSIDERATION  OF  INITIAL  IMPERFECTIONS 

II 1. 1,  Introduction 

The  Euler-Bernoulli  beam  type  theorem  assumed  in  section  II. 1  which  is 
similar  to  the  Kirchhoff  assumption  of  plate  theory  or  Love's  first  approxi¬ 
mation  of  shell  theory  results  in  neglecting  the  transverse  shear  deformation 
and  transverse  normal  stress. 

Due  to  recent  developments  of  new  types  of  material,  the  orthotropic 
feature  of  many  materials  must  be  taken  into  account.  The  present  investi¬ 
gation  is  particularly  concerned  with  materials  in  which  the  reinforcement 
consists  of  glass  fibers  embedded  in  an  epoxy  matrix  with  every  glass  fiber 
oriented  toward  the  center  of  curvature  of  the  structure,  in  this  case  an 
arch.  This  radial  filament  construction  offers  the  advantage  of  the  fiber 
being  loaded  in  transverse  compression  rather  than  in  axial  compression 
which  may  lead  to  column-type  instability.  The  resin  is  restrained  from 
extrusion  from  between  the  glass  fibers  by  shear  coupling  to  the  fibers, 
thereby  placing  the  filaments  under  axial  tension.  Thus  the  inclusion  of 
transverse  normal  stress  and  shear  effects  may  be  of  importance.  With  inclu¬ 
sion  of  these  effects,  an  Euler-Bernoulli  beam  type  assumption  is  inade¬ 
quate.  And  a  modification  to  Euler-Bernoulli  type  assumption  (c)  lies  in 
not  completely  neglecting  the  ratio  z/R  with  respect  to  unity  but  instead 
writing  j+zJr  '  1  -  £  *  %)  ™  the  strain-displacement  relations. 

The  initial  radial  imperfections  of  the  middle  surface  are  included  in 
the  derivation  to  more  nearly  agree  with  reality.  This  leads  to  a  lower 
asymmetric  buckling  load.  Reissner's  Variational  Theorem  [21]  is  used  to 
derive  the  equilibrium  equations  as  well  as  the  stress-strain  relations. 
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The  numerical  method  outlined  in  the  foregoing  chapter  is  used  to  solve  the 
equations 


III. 2.  Strain-Displacement  Relations 

Let  us  use  the  coordinate  system  described  in  section  II. 2.  and  also  the 
von  K&rman  type  non-linear  strain-displacement  relations.  The  transverse 
normal  strain  and  transverse  shear  strain  yv,  together  with  e  now  des- 

E  aZ  X 

cribe  the  state  of  strain,  viz: 


:  =  — rr-  (U,  +  jr)  +  - (^  W»2  +  W,  W  y) 

x  i+f  x  R  (i+f)2  2  x  x  °*x 


cz  =  W*z 


(3.1) 


Yx,  =  -rr  W  +  U,„ - -  u 

XZ  1+|  *  R{1+|) 

where  wQ  is  the  initial  imperfection  in  the  radial  direction  (see  [18]).  Let 
us  introduce  the  displacement  approximations 


U  (x,z)  =  u  (x)  +  ?v'  (z) 

W  (x,z)  =  w  (x)  +  zv*  ‘  (x)  +  ^  z2  w"  (x) 


(3.2a) 

(3.2b) 


Substituting  (3.2)  into  (3.1)  and  multiplying  (3.1a)  and  (3.1c)  by  (l+j|),  we 


(,+f)  ex  *  Ex *  zKx +  f  +  fr  K  *  fr  Kx' +  It  kxv  +  fr  Kx 


R'  ‘X  ‘X  •  X  2  ’ X 


*z  ■  +  ZKZ 


(3.3) 


<1+t)  Txz  =  Txz  *  zKxz  +  r  Kxz 


where 


=  U.  +  F  W  +  i  W,2  +  W,v  W. 

X  A  K  Z  A  X  0  }  A 

i  1  •  Ci  2  1  •  C] 

Kx  u,x  +  R  w  2R~  w’x  +  w’x  w*x  +  w,x  wo,x  "  R~  w*x  wo,x 


(3.3d) 


(3.3e) 


,  ,  ,2  „  2C,  •  ci  2  2C] 

u  =  F  W"  +  W,v  +  W,v  W,v  -  -5-  W,v  W,v  +  T  W,  +  -5-  W»  W. 


'X  R 


■x  w’x  '  —  w’x  w*x  +  7  w*x  T  JT  w»x  "o»x 


(3-jf) 


TT  wo,x  w’x  +  wo,x  w*x 


11  -I  1  11  C-.  1 2  -i  ti  C-.  1  C-j  1 

Kx  =  6>(2w*x  w’x  "  2R  w*x  ”  2R  w*x  w’x  *  ^2  w*x  w’x  +  ^2  wo,x  w*x 

C-j  11  (3.3g) 

”  2r  wo,x  w>x> 

n  C-j  I  11  C-j  1 2  C-j  It  C-j  II 

k"'  =24-(1w,^  "  “w*x  w,x  +  ^2  w»x  +  ^2  w*x  w*x  +  ^2  wo,x  w,x*  (3*3h' 


iv  C-j  11  C-j  1  11 

Kx  =  120‘(“§Rw*x  +  ”2“w*xw'x) 

K  =  720  *( — x-  w,  ) 
x  8R^  x 


0 

ez  =  w 


(3.3i) 


(3.3 j) 


(3.3k) 


Kz  =  w 

^xz  =  -t+u+w-x 


Kxz  =  w’x 

I  M 

and  Kxz  =  w,x. 


(3.31) 

(3.3m) 

(3.3n) 

(3.3o) 


The  coefficient  C-,  is  introduced  to  distinguish  the  modification  in  this  theory 


from  the  Euler-Bernoulii  type  theorem. 


III. 3.  Derivation  of  Stress  Strain  Relations  and  Equilibrium 
Equations 

If  the  force  resultant  N,  the  moment  M  and  the  transverse  shear  resultant 
V  (all  per  unit  length  of  middle  surface)  are  defined  in  the  conventional  manner, 


h  h  h 

'2  ( 2  [Z 

N  =  o  dz  ,  V  =  t  _  dz  ,  M  =  o  zdz. 

h  x  h  xz  J_h  x 

~Z  "2  2 


(3.4) 
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then  we  az>z>ume  that  (see  [13]) 


_  N  ,  H  2  _  N  12M  _ 

°x"h  h2/6h/2-h  ~^z 

.  3  Vr,  ,  z  >2-  3V  r.  4z2! 
Txz  "  ?  IP-1  ~  V2'  J  "  2h  L1  "  h2J 


(3.5) 


(3.6) 


*1+f)  °z  =  C2  ^2  h  +  4h  ^h/2^1  "  ^h/2^  +  2  ^^h/Z1  "  Vhjzf^  *3*7* 

Here  the  functions  S  and  T  are  as  yet  undetermined;  q  is  the  value  of  the 
load  applied  at  the  top  surface  of  the  arch, 

and  C2  is  introduced  for  the  purpose  of  distinguishing  between  the  contribu¬ 
tions  of  transverse  shear  deformation  and  normal  stress,  (C2  =  0),  otherwise, 
it  should  be  regarded  as  unity. 

A  convenient  derivation  of  appropriate  stress-strain  relations,  equili¬ 
brium  equations  and  natural  boundary  conditions  can  be  carried  out  by  means 
of  a  variational  theorem  due  to  Reissner  [21],  viz: 


6  {///  [  \  Ex  +  Txy  Yxy  +  “z  Ez  -  A]  (,Ht)  dx  dy  dz 

2  C  h 

-  //  q(w  +  ~  w  +|  w  )  (1+  ^-)  dx  dy}  =0 


(3.8) 


liittei  w  •  VI  U||  VI  VltUbl  vpl\.  lllU  UCI  IUI 


,  1  r  2  .  *xz  2  2x  2  9  •, 

A  -  2E^  Lox  +  ^  °z  +  r  T«z  -  2vxz  °x  °z]- 


(3.9) 


Integrating  (3.8)  with  respect  to  z,  carrying  out  the  variations  and 
integrating  by  parts,  the  stress-strain  relations,  the  equilibrium  equations 
and  associated  boundary  conditions  are  obtained. 

t  If 

The  coefficients  of  and  <5w  give  the  expressions  for  S  and  T  res¬ 
pectively.  A  simplification  is  made  by  neglecting  the  nonlinear  terms  and 
also  the  shear  term  V,  since  as  the  arch  becomes  thinner  the  effects  of  S  and 


o  * 


T,  i.e.  the  effects  of  transverse  normal  stress,  becomes  negligible.  Thus 
the  equation  (3.5)  In  reference  [13]  is  used: 


S  =  .  M 

5  R 


T  -  5h  M  ,  3hH  n 

T--opN+-5-q 


(3.10) 

(3.11) 


The  coefficients  of  sll  and  6M  give  the  relations  for  M,  N  and  displace¬ 
ment  components.  They  are 


N  ■  Exh  (AA)  [1  -  (1  ♦  t  Exh  •  CM  •  (BB) 

.  _L|L_  [caa  .  |  (1  +  -  3  •  CM  +  h] 


(3.12) 


f:„b  C„v 
f  ft  j.  2  xz 


M  =  CA  •  {-—  (1  +  -^)  (AA)  -  Exh  •  BB  +  Hq  [|  (1  +  ^^i)  -3]}  (3.13) 


where 

CA  =  l/f^l  +  ^JS)  (1  ♦  C2vxz)  -  If] 
CAA  =  CA  •  (1  +  C2vxz)/R 

-  c°  +  u*  .  h4  «»  ,  h6  „v 
M  Gx  T2  Kx  4J-80  Kx  6T'4"48  Kx 

rr  =  v  j.  _3h^  i,"  ,  3h^  ./iv 

Kx  3!-20Kx  +  5Tn2Kx 

The  coefficient  of  <5V  gives 

V  .  ^  (v?.  ♦  t  K.’.  J  -  m  •  CC 

«  AA  tU  A£  O 


1A\ 
• ' / 


where 


o  1 

cc  ■  4  +  4o 


Equations  (3. 1C)  through  (3.14)  are  the  stress  resultant-displacement 
relations. 

The  coefficients  of  ou,  6u  ,  6w,  6S  and  6T,  in  this  case  yield  the  five 
equilibrium  equations.  These  are 


N,x  +  R  =  0 


(3.15) 


M,x  -  V  =  0 


(3.16) 
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1  1.2  it  h  C-j  i  C,h2  C^h  ii 

N  C"  R  +  w’xx  +  wo,xx  +  2?  w’xx  “  12R  w‘xx  +  ^  (w*xx  +  wo,xx)  +  :^2w>xx] 

h2  „  C,h2  C,h2  ,  C-, h4  „  ' 

+  N.xCw.x  +  wQ>x  +  24  w,x  -  Tg F  w»x  +  ^2  (W*X  +  Wo,XJ  +  ^2  ^ 

C,  ,  3C,h2  „  3C,h2  , 

+  M  (w,xx  +  wQjXX)  +  w»xx  "  4or  w»xx  +  40r2""  w*Xx] 

C,  ,  3C,h2  „  3C,h2  , 

+  M*x  R~  ^w*x  +  wo,x^  +  w*x  “  TOfTw*x  +  4oi?~W,x'* 


+  V.x  +  q’H  =  0 


(3.17) 


(r+^V-)  - 


JL  c  u  n  h2  \  +  H  C2vxzh  /6  ^2 

E  h  ^2vxz  *'  5  '  +  Ru  '7 

txn  ^  XZ  28R  v  E  h2  Rvzx  5  70R' 

fcA  A 

h2C, 


M  -2vxz"  ,6  3h"  x  C2vxz  ,,  u,,ul 


^7zx 


6hC, 

o-asr 


+ - Hq  +  c  e°  =  G 

20R2  2  2 


(3.18) 


N  C2vxzh  ,  1  ,  h2  %  ,  H  C2vxz  n 
Exhl2Ryzx(?  84R2)  i?  5  0 

A 


14R2v 


-)  - 


l2  Cnh 

+^?]Kq+^=° 


zx 


C2vxz  r4  h_ 
20Exvzx  L?  -  6R 


(3.19) 


From  (3.16) 


V  =  M,, 


(3.16a) 


Substituting  (3.16a)  into  (3.15)„  we  get 


“•x  *  if  H'x  ■  0 


(3.15a) 


or 


N,x  =  "  R  M’x 


(3.15b) 


Substituting  (3.15b)  and  (3.16a)  into  (3.17),  we  get 

H  *  C-  7  +  (CD) » x3  +  M,x  •  [(CE)  -  1  (CD)]  +  M  •  [(CE),X3  +  M,xy  +qH  =  0 

(3.17a) 


whei  e 


u2  n  h  Ci  i  C-«h2  C-ih^  ii 

CD  =  w'x  +  wo.x  +  2?w-x  -  Ww'x  +  ^  <w'x  +  wo,x>  +  ^2  w-x 
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C,  ,  3C,h2  „  3C,h2  , 

CE  =  -  r  <*s +  wo,x> *  w>x  ■  tot  w>x +  w*x 

Equations  (3.15b),  (3.16),  (3.17a),  (3.18)  and  (3.19)  are  the  equili¬ 
brium  equations. 

The  natural  boundary  condition  terms  are 


[[6U  [N]  +  6u'  [M]  +  6w  [(N  +  — \  N  -  l  M)  (w,x  +  w 0>x) 

1 2R  * 


C,h2  3C,h2  ,  ,2  3C,h2  C,h4 

+  (H  “  W-N  f-lj-M)  w,x  +  -Tfe-M  +  ~^N)  W,x] 

,  .2  3h2C,  .2  -,.2  .4  , 

+  6w  [(M  -  jgl  +  M)  (w,x  +  wQ>x)  +  (^4  N  "  20"  M  -  §0R  N)  w>x 


ol.2  C,h4  3C,h4  „  „  ,2  3C,h2  C,h4 

+  <frH  “’x]  +  5W  [(27n  •  Tor"  h  +  ^2n) 

,,2  C,h4  3C,h4  ,  h4  „  3C,h4 

^W,X  +  wo,x^  +  * W  M  ■  T60R  N  +  ^2  w,x  +  320  N  W,x  “  338R~  M 


C-jh  „  h2  , 

+  -J-T  V  w>x  +  W  V]  ]]o 

3584R^  X  w  0 


(3.20) 


They  are  in  the  form  of  stress,  displacement  or  mixed  boundary  conditions. 
For  example,  we  have  for  the  first  term,  either  u  prescribed  or  N  =  N  (constant) 
prescribed  at  x  =  0  and  L. 


III. 4.  Non-dimensional i zed  Equations 

Let  us  introduce  the  following  non-dimensiona1  variables: 


x  =  Is, 


h2 

u  =  r-  4. 


W  =  ha, 


«'  =u6  > 


u"  -  h 
W  -  —9y 

l6 


w  =  h  wn 
0  0 


x  =  L7Rh 


and  (jj-)n  =  pn  where  n  is  en  integer. 


Then  the  abbreviated  constants  and  functions  are 


I 


CA  =  —  = _ _ _ _ 1  _ 

h  pV-o  +  -£“)  (i  +  c2„sz) .  12 


A-pg-(1  t  C2vxz) 

p4,»2-o  +  -^)  (i  +  c9v  j  - 


H  =  1  +  ^  •  A  .  p2 


(3.21) 


(3.22) 


(3.23) 


=  7  =  V  +  A  •  «  +  r  o.f  +  w0lS  f.s  +  ar  P2  •  (Ay  +  s,2 
+  «,s  1,s  -  2C,  A  •  P1  •  a,s  B,s  +  C,  A2  •  p2  a,2  +  2C,  A2  •  p2 


'  Vs  °>s  '  2C1  '  A  •  p’  •  W0>s  8,s  +  «0>s  Y.s)  +  gj  p4  •  y.j  -  r  1 


1  ^1  9  2  ?  9  9  C,  o 

*  P  ‘  B»s  Y»s  ~  21  X  P  ‘  P*s  +  F*  ^  *  pZ  *  a»*  Y>~  +  ^  X2 


*s  Y,s  T 


'  p2  ■  Vs  Y’s>  +  m  •  P8  •  V.? 


(3.24) 


dr-  _  h-BB  .  ,  .  ,  1  ci 

BB  ~  —  -  ##s  +  a  •  p  *  9  -r 


2  2.1 


2  ^>s  +  A  •  p  0  -  tj—  A  •  p  •  a,s  +  p  •  a,s  0,s  +  p 


•  Vs  8*S  -  C1  •  1  •  p2  •  Vs  “’S  +  ^  •  p3  •  <?  6>s  T.s  -  r  x 

1  2  1  o  o  o 

P  *  3>_  -  A  •  p  •  a,.  Yi.  +  C,  k  •  X)  •  ft.  +  r.  \2 


2”  A  •  p  ’  a’s  Y’s  +  C]  X  *  p  *  a*s  0*S  +  C1  x 


*  P  *  w0,s  6’s  '  2]  x  *  P1  *  w0,s  Y’s>  +  TT2  p5  ■  (“  r*  Xy*s 
^1  ? 

*  2~x  •  *»s  Y’s} 

CC=^=-p2.  A  •«+  *  4-  B,s  +^P2  Y,s 

CD  =  7 =  a»s  +  wo,s  +  2? p2  *  Y’S  -  re  * x  •  p3  •  e»s  +  4 x2  *  p4 


'1  s.  2 


(3.25) 


(3.26) 


'  K  +  wo,s}  +  160  ^  *  P6  *  Y’< 


(3.27) 


fir  _  CEJi  _  1 

Lt  y—  -  p 

P 

3C,  o 


1  =  p  C"  C1  x  pl  '  (a»S  +  wo,s)  +  6»s  "  W~  x  *  P3  *  Y’s 


+  ZL  >2  4  _ 

+  40  x  '  P  *  8>s] 


(3.28) 
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and  the  stress  resultants  are 


G 

N  8  IT  =  m  '  P?*  •  V  -  x  *  CAA  •  p2  .  (1  +  -OZ)]  +  COT  •  BF  •  p2 
-  •  [COT  *  p2  *  X  •  (1  +  ~^)  -  3  •  COT  +  1]  #. 


=  BF  •  p2  •  OT  +  GOT  •  p2  •  BF  -  OF  •  q 


M  =  —^2  =  cff  {p4  •  X  •  (1  +  AA  -  p2  •  BB  +  —p.  H 

Eh  6  z 

A  Q 

J  [p2  •  X  •  (1  +  -  33  /-} 


=  DC  •  p2  •  AA  -  CF  •  p2  *  BB  -  DE 


V  =  OT  =  6E  pl  *  p2  '  x  *  *  +  *  +  °*s  +  40  Y*s] 

A  X 

=  DF  •  p1  *  CC 

—  S  ?  _ 

S  -  c  L  =  “  p  X  •  M 
X 

T“rh  =  -f  P2  *  A  '  «  +|h  •  q 


(3.29) 


(3.30) 


(3.31) 


(3.32) 


(3.33) 


where 


O 

DA  =  1  -x  •  CAA  •  p2  •  (1  +  ^r5-) 

DF  =  {COT  •  [p2  -  X  .  (1  +  ~^)  -  3]  +  1}  •  F 

£  *  C 

DC  =  Cff  •  p2-  A  ■  (1  +^Si) 

_  _  CflV  n  C„V 

DE  =  -  CA  •  H  •  [p2  •  X  •  (1  +  -  3] 


a  = 


DF  =  ~ 


With  these  expressions,  the  non-dimensional  equations  are: 


9^^ 


mm 


m 


(i)  DG  •  (M),  +  OFT  •  (bB),  =  0 

«>  s 

whe.’e 

DG  =  DA  +  X  •  p2  *  DC  and  DiT  =  CKK  -  p2  •  x 

(ii)  t)C  •  p2  •  (M),s  -  Cft  •  p2  •  (BB),S  -  DF  •  CC  =  0 

(ill)  SA  •  [OS'  *  (-  X+  CD,$)  +  DC  *  CE,$]  +  BB  •  [C®-  •  (-  X  +  £d,s) 
-  W  ‘  ^*$]  +  M,s  •  DC  •  (CE  -  X  •  p2  •  W)  -  BB,S  •  ZK 

•  (CE  -  X  •  p2  •  CD)  r  M",ss  •  DC  -  BF,ss  *  CK  +  [H  -  DB 

•  p2  •  (-  X  +  CD ,  )  -  CF  •  p2  *  CB, J  -0=0 
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(3.34) 


(3.35) 


(iv) 


vv; 


where  Q«-| 

DL*M'+DM*BB  +  DN-Q+0  =  O 
whe^e 

DL  =  f  P~  •  DA  +  DJ  •  DC)  •  p1 

DM  =  (15T  •  CKK  -  DJ  •  CK)  •  p^ 

DN  =  (DK  -  DB  •  DT  -  DJ  •  DE)  •  p 
_  C?x2  . 

D!  ■  C2  '  wxz  -  ~w  P ■  PR  „ 
—  3  ?  -4 

ni  ^  J  /»  -n 


°J  =  |  C2  ■  pR  •  pZ  •  A  •  (2  +  ■  A2) 

6x 

35  K  "20 


—  up  — 
DK  =  tj=-  H 


pR  •  (-  1  +  |i  •  pZ  -  i2 


pR  =  vxz/',zx 


ftft 

LO 


TP  .  pr=-  rrr  _2 

«H  T  Ul  *  BB  -  UU  *  H  *Q+Y=0 


where 


^  ■  !2cr  •  vxz  •  <’  -  r£~  »2> 

zx 

DR  =  -/pR  •  H  (-  f  +  A  •  p2  -  A. 

DS"  =  DP"  •  DA  +  DQ"  •  CA 
DT  =  DP"  *  CM"  -  DQ*  *  CTT 
'DU  =  DP  •  DB  +  DQ  ■  DE  -  DR 


p4) 


DP  =  Cg  •  pR  •  p  •  X  •  (f-  +  -gg  p^  •  x2) 


(3.35) 


(3.37) 


(3.38) 
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In  deriving  the  last  two  equations  (3.37)  ana  (3.38),  we  take  C 2  -  1  (tor 
the  last  terms  in  both  equations)  in  order  that  the  equations  are  nonvanish¬ 
ing  (i.e.  the  coefficient  matrix  is  nonsingular).  This  is  for  the  case  of 
neglecting  transve. >e  normal  stress  (C2  =  0)  and  the  solutions  $  =  0  and 
Y  =  0  are  obtained. 


III. 5.  Comparison  of  Three  Theories 

When  we  set  C2  =  6  =  y  =  0  in  (3.34)  through  (3.38),  the  equations 
are  specialized  into  those  of  shear  deformation  theory  with  the  modification 
in  the  Euler-Bernoulli  type  theorem.  If  we  completely  neglect  z/R  terms 
(C-j  =  0)  as  well  as  imperfection  terms  (wQ  =  0)  the  following  shear  theory 
is  obtained? 
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(0  -  (A  -  «,ss)  (♦,,  +  A  ■  a  +  £  «,f)  -  ^  p2  '  A  •  a>s  •  *.ss 

+  TT  **sss  -  Q  =  0 

The  underlined  terms  are  the  differences  in  these  two  sets  of 

equations.  However  the  underlined  terms  all  contain  powers  of  h/L, 

2  4 

i.e.,  p  ,  p  ,  etc.  Thus  the  differences  between  the  two  sets  becomes 
smaller  as  the  arch  becomes  thinner.  Further  simplification  of  the 
equations  as  h/L  approaches  zero,  leads  to  the  classical  equations 
for  the  arch,  viz: 

(a)  4,ss  *  A  •  o,s  +  a>s  a,ss  -  0  (3.42) 

(b)  ^  +  a>s  =  0  (3.43) 

(c)  (*,s  +  A  a  +  1“,^)  (A  -  a,ss)  -  ^>sss  -  Q  =  0  (3.44) 

With  the  initial  imperfections  neg1  "ted,  the  first  equation  is  identical  to 
(2.18a)  while  combination  of  the  last  two  equations  is  identical  to  (2.18b). 


1 1 1. 6.  Method  of  Solution 

Let  us  employ  the  same  method  of  linearization  of  the  nonlinear 
equations  as  in  the  previous  chapter.  For  boundary  condition,  the  clamped 
end  case  of$  =  ^  =  a=  6  =  Y  =  0ats  =  Q  and  1  are  considered,  since  these 
displacement  components  can  be  rather  easily  specified.  Because  third  order 
derivatives  are  involved  in  the  third  equation,  forward  and  backward  diffe¬ 
rences  are  used  in  order  to  avoid  the  values  at  nodal  points  outside  the 
boundary,  viz: 


43a. 
S  1 


_  D' 


("  3cti. 


+  10a.. 


12ai+i 


*  6ai+2 


"  cti+3) 


and 


7sS 


2 'ai-3 


6ai-2  +  12ct  i-1 


lOa^  +  3a^+1) 
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where  A  is  the  forward  difference  and  v  is  the  backward  difference. 

As  in  the  last  chapier,  the  equations  are  written  in  finite  difference 
font?.  Deflection  input  in  the  computer  program  will  be  used  to  extend  the 
load-deflection  curve  beyond  the  first  peak.  Therefore  the  additional  equa¬ 
tion 

AREA  =  l  a. 

N+l  i=l  1 

is  used. 

The  newly  formed  equation  represents  a  system  of  (5N  +  1)  simultaneous 
nonlinear  algebraic  equations  which  are  solved  by  Newton's  iteration  scheme. 
The  polynomial  in  section  11,4.  is  also  used  as  for  the  antisymmetric  initial 
imperfection. 

The  orthctropic  elastic  constants  are  taken  from  reference  [22]  as 
Ex  =  4.20  x  106  psi 
Ez  =  10.20  x  106  psi 
G  =  2.7  x  106  psi 

vxz  “  O*11 
vzx  -  °*26 

where  the  relation  E  v  =  E  v  is  approximately  satisfied. 

A  4a  4  Am 

1 1 1. 7.  Discussion  of  Results 

Theories  of  varying  degrees  of  accuracy  may  be  obtained  as  special 
cases  of  the  general  theory  which  considers  transverse  normal  stress  and  shear 
deformation.  For  the  limiting  point  "a"  as  well  as  the  bifurcation  point  "b" 
(Fig.  8),  the  more  accurate  the  theory,  the  lower  the  critical  load  is  found 
to  be.  The  load-central  deflection  curves  for  \  =  56.7  are  shown  in  Fig.  8 
and  their  enlarged  upper  portions  are  shown  in  Fig.  9.  The  discrepancy  in 
critical  load  for  the  classical  theory  developed  in  the  last  chapter  and  the 
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(in  CHAPTER  II) 

CLASSICAL  (h/L*0) 


irr  a  r>  - 

onc«n  ucrunivj.  i 

(in  APPENDIX) 
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one  obtained  as  a  special  case  from  the  general  theory  is  about  1.3  percent 
of  the  latter.  This  discrepancy  is  a  result  of  a  different  network  of  grid 
points  involved  in  the  finite  difference  formulation  of  the  fourth  order 
derivative  in  equation  (2.18b)  and  the  combination  of  first  a'.id  third  order 
derivatives  in  (3.43)  and  (3.44)  respectively.  For  comparative  purposes, 
all  theories  will  be  compared  to  the  classical  type  theory  obtained  as  a 
special  case  from  the  general  theory  which  includes  transverse  shear  as  well 
as  transverse  normal  stress  effects.  For  h/L  =  1/6G,  the  asymmetric  critical 
load  corresponding  to  classical  theory  is  about  4.3  percent  higher  than  that 
found  when  the  transverse  normal  stress  and  shear  deformation  theory  is 
employed.  Also  the  load  found  from  modified  shear  deformation  theory 
(C2  -  8  =  y  =  0)  is  only  0.7  percent  greater  than  the  general  theory.  The 
shear  deformation  theory  (C-j  =  C2  =  B  =  v  -  0)  runs  1  percent  higher  than 
the  modified  shear  deformation  theory.  Similar  results  are  obtained  for 
A  =  35.0  as  shown  in  Fig.  10. 

The  shear  effects  become  very  significant  when  the  arch  is  moderately 
thick.  For  X  =  13,  h/L  =  1/13.3  as  shown  in  Fig.  11,  the  critical  load  found 
from  classical  theory  is  16  percent  greater  than  that  found  when  normal  stress 
and  shear  are  considered.  However  *ne  difference  between  results  obtained 
from  the  modified  shea*'  theory  end  the  normal  stress  and  shear  theory  is 
usually  less  than  2  percent.  In  this  case  the  shear  deformation  theory  (with 
Ci  =  C2  =  8  =  y  =  0)  is  not  adequate  to  describe  the  behavior  of  the  arch 
since  the  critical  load  is  about  8  percent  higher  than  that  found  from  the 
modified  shear  deformation  theory  (C-j  /  0). 

The  comparison  of  the  critical  loads  found  from  the  classical  as  well  as 
the  shear  and  normal  stress  theories  is  shown  in  Fig.  12.  The  percentage 
difference  between  results  obtained  from  the  former  and  the  latter  theories 


is  noted. 


\  *  \ 
\\\  \ 

V  \ 

\\  \ 

\\ 

\\\ 


■ - CLASSICAL  (  h/L=0) 

- SHEAR  DEFORM., Ci=0 

( h/L=  1/50) 

• - SHEAR  DEFORM.,  | 

(h/L=  1/50) 

— -  GENERAL,  Cj=C2- 1 
(h/L=  1/50) 


0.9  1.0 

-wc/h 
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CHAPTER  IV 

EXPERIMENTS  ON  RADIALLY  REINFORCED  ARCHES 
IV.l.  Introduction 

Experiments  v/ere  carried  out  on  arches  composed  c~  resin  matrices  in 
which  were  embedded  small  diameter  glass  fibers  each  oriented  toward  the 
center  of  curvature  of  the  arch.  This  type  of  construction  was  first  pro¬ 
duced  for  laboratory  purposes  by  Uniroyal,  Inc.  [23],  The  nature  of  this 
composite  material  makes  it  well-suited  to  withstand  loadings  normal  to 
its  convex  surface.  Such  loadings  would  occur  in  deep  submergence  pressure 
hulls.  Evidently  during  hydrostatic  loading  both  the  resin  matrix  as  well 
as  the  glass  fibers  are  acting  in  compression  and  the  resin  is  restrained 
from  extrusion  from  between  the  glass  rods  by  shear  coupling,  thus  placing 
the  glass  filaments  in  a  state  of  axial  tension.  It  has  been  demonstrated 
that  spherical  shells  constructed  of  such  radially  oriented  reinforcement 
nave  extremely  desirable  weight/displacement  relations  [24]. 

The  objective  of  the  present  experiments  is  to  investigate  validity  of 
the  various  theories  developed  earlier  in  the  present  investigation. 

IV. 2.  Fabrication  ef  the  Material 

Fabrication  of  the  radially  .-einforced  arches  was  carried  out  in 
accordance  with  the  following  procedure: 

(a)  Material:  The  material  utilized  high  strength  S  glass  roving 
embedded  in  an  epoxy  resin  which  is  marketed  through  Minnesota 
Mining  and  Manufacturing  Co.  under  the  designation  of  1009-26-S 
which  is  a  glass  impregnated  tape  one  inch  v/ide  and  SCO  feet  in 
length.  The  glass  fibers  are  of  course  oriented  in  the  500  foot 
direction.  The  tape  is  cut  into  short  lengths  each  corresponding 
to  the  desired  thickness  of  the  arch,  in  this  case  0.25  inches. 
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(b)  Curing:  To  build  sr.~11  mosaic-type  curved  blocks,  the  0.25  inch 
lengths  of  tape  were  stacked  into  a  curved  mold  having  a  16  inch 
radius  of  curvature.  A  spring  force  in  the  direction  of  the  arcn 
middle-surface  was  applied  so  as  to  develop  a  pressure  of 

2000  lb/in^  to  compact  the  short  lengths  of  tape.  Holes  were 
provided  in  the  mold  to  permit  escape  of  excess  resin.  The  mold 
was  then  placed  in  a  small  oven  and  the  mosaic  block  cured  at 
temperatures  of  300°F  to  325°F  for  30  minutes.  The  cured  block 
contained  about  80  percent  glass. 

(c)  Post-curing:  The  cured  mosaic- type  blocks  v/ere  joined  to  one 
another  by  an  epoxy  cf  minimum  thickness  so  as  not  to  influence 
the  resin  content  of  the  resulting  arch.  The  entire  assembly 
was  tnen  placed  in  a  larger  mold  and  cured  again  at  350° F  to 
3/5°F  for  four  hours. 

In  addition  to  arches  constructed  by  the  above  technique,  additional 
specimens  were  provided  by  Uniroyal,  Inc.  The  physical  properties  of  the 
resulting  composite  are  described  in  [22]  and  Fig.  13(a)  illustrates  the 
nature  of  the  material. 

IV. 3.  Description  of  Test  Apparatus 

The  test  apparatus  consists  of  a  loading  fixture  which  includes  end 
supports  for  the  arch  as  well  as  a  system  to  provide  lateral  pressure  on 
the  convex  surface  of  the  arch.  Also,  a  system  of  transducers  for  measuring 
displacements  was  developed. 

(a)  Support  conditions:  The  ends  of  the  arch  were  clamped  by  means  of 
two  metal  blocks  tightened  around  the  arch  by  means  of  two  bolts. 
One  of  these  bolts  went  through  the  end  of  the  arch.  It  was 
believed  that  this  corresponds  reasonably  well  to  a  clamped  end 
condition,  as  illustrated  in  Fig.  13(a).  A  slot  corresponding  to 
the  curvature  of  the  arch  was  milled  into  the  blocks  so  as  to 
provide  good  clamping.  This  technique  for  simulating  a  clamped 
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end  was  first  tested  for  a  straight  cantilever  beam  and  found 
to  lead  to  deflections  in  good  agreement  with  simple  beam  theory. 

On  this  basis  it  was  considered  a  good  representation  of  a  clamped 
end  condition. 

(b)  Deformation  measurement  system:  Radial  displacements  were  moni¬ 
tored  by  use  of  three  Schaevitz  Linear  Variable  differential  trans¬ 
formers,  Type  300-SSLT,  henceforth  designated  LVDT.  These  units 
each  has  a  range  of  t  0.3  inches  with  a  linearity  larger  than  i  1 

percent  of  the  full  range  output.  Each  LVDT  consists  of  a  three 

winding  transformer  together  with  a  removable  core  as  illustrated 
in  Fig.  14(a).  When  the  core  is  centered  with  respect  to  the  two 
secondary  windings  the  induced  voltages  V-j  and  are  equal.  As 
the  core  is  displaced  toward  either  secondary  winding  that  induced 
voltage  increases  whereas  that  from  the  other  secondary  winding 
decreases.  When  the  core  is  attached  to  the  arch  the  radial  dis¬ 
placement  of  the  arch  can  be  monitored.  With  the  aid  of  some 

additional  circuitry  the  interpretation  of  resulting  data  can  be 
made  easier.  The  circuit  corresponding  to  this  is  shown  in  Fig. 
14(b).  With  the  core  in  the  central  position  the  output  voltage 
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increases  positively  and  vice  versa.  Fcr  full  core  displacement 
of  t  0.3  inches  the  output  voltage  changes  by  approximately  t  3 
volts  for  a  6.3  volt  excitation  of  the  primary. 

The  arch  displacements  treated  were  of  the  order  of  magnitude 
of  the  arch  thickness,  i.e.  0.25  inches.  It  was  anticipated  that 
the  maximum  radial  displacement  would  be  no  more  than  0.5  inches. 
Because  of  this  the  core  was  attached  to  the  arch  after  a  positive 
displacement  of  G.25  inches  of  the  core  had  been  imposed.  This 
leads  to  a  useful  dis'lacement  range  of  0.50  inches  whicn  is  within 
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(b)  ADDITIONAL  CIRCUITRY  TO  MAKE  INTERPRETATION  EASIER 

F,ai4  LINEAR  VARIABLE  DIFFERENTIAL  TRANSFORMER  a  CIRCUITRY 


the  linear  range  of  action  of  the  LVDT. 

The  output  voltage  vQ  was  monitored  by  a  Honeywell  Visicorder, 
Type  2106  which  displayed  voltage  changes  on  a  six  inch  wide  strip 
chart.  The  signal  conditioning  amplifiers  were  adjusted  >o  that 
a  core  displacement  of  0.5  inches  (measured  from  the  pre-set  0.25 
inch  position)  gave  a  full  five  inch  displacement  on  the  chart 
paper.  Thus,  0.1  inch  of  deflection  on  the  chart  corresponds  to 
0.01  inches  of  arch  deflection. 

(c)  Pressure  loading  system:  The  pressure  system  consists  of  an  oil 
jack,  a  pressure  chanter,  and  a  pressure  gage. 

(1)  Pressure  chamber:  This  was  formed  on  the  convex  side  of  the 
arch  by  bounding  that  side  by  a  concave  shaped  piece  of  hard 
wood  together  with  two  0.5  inch  thick  plexiglass  sheets,  one 
on  either  side  of  the  arch.  The  plexiglass  sheets  were  bolted 
together  through  the  hard  wood  as  well  as  through  the  arch 
supports.  A  bicycle- type  inner  tube  was  placed  in  the  space 
between  the  arch  and  the  hard  wood  and  connected  to  the  oil 
jack.  Uniform  hydrostatic  pressure  was  then  applied  to  the 
convex  side  of  tne  arch  by  pumping  oil  into  the  extremely 
flexible  inner  tube. 

(2)  Pressure  gage:  A  gage  manufactured  by  the  U.S.  Gage  Co.  and 
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capable  of  monitoring  pressures  up  to  200  Ib/in  was  connected 
between  the  inner  tube  and  the  oil  jack.  The  gage  was  of 
course  calibrated  prior  to  testing  the  arches. 

(3)  Oil  jack:  A  Hein-Wemer  Pushmaster  Hydraulic  Pump  Unit  Model 
FP-4  was  employed.  It  was  found  that  ordinary  jack  oil  rather 
quickly  damages  the  rubber  in  the  inner  tube  so  paraffin  oil 
was  employed. 
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IV .4.  Test  Procedure 

The  cured  arch  was  clamped  in  the  metal  blocks  and  this  assembly 
surrounded  by  the  plexiglass  plates  and  hard  wood  backing  for  the  Inner 
tube.  The  core  of  each  LVDT  was  epoxied  to  the  concave  surface  of  the  arch. 
One  transducer  was  located  at  each  quarter  point  of  the  arch  and  one  at  the 
mid-point.  Figs.  13(b)  and  (c)  illustrate  the  test  apparatus. 

The  output  sides  of  the  LVDT's  were  connected  to  the  input  side  of  the 
Honeywell  Accudata  117  DC  Amplifier  and  the  LVDT's  energized  by  means  of  a 
transformer  operating  from  the  115  volt  AC  source.  The  amplifier  was  acti¬ 
vated  and  the  ultraviolet  light  beam  centered  on  the  chart  paper.  The  gain 
was  adjusted  so  as  to  keep  the  trace  on  the  chart.  With  the  aid  of  a  micro¬ 
meter  the  LVDT  was  displaced  in  the  direction  of  the  convex  side  of  the  arch 
by  0.25  inches  and  the  LVDT  locked  in  position. 

Oil  pressure  was  slowly  applied  while  the  Visicorder  chart  paper  was 
driven  at  0.5  inches/sec.  LVDT  readings  were  taken  at  10  lb/in  increments 
until  arch  collapse  occurred. 

IV. 5.  Test  Results 

Eleven  arches  were  tested  and  the  results  are  indicated  in  Table  1 
and  plotted  as  load-deflection  diagrams  in  Fig.  15.  This  permits  compari¬ 
son  between  the  various  theoretical  results  obtained  earlier  in  this  study 
and  experimental  evidence. 

Specimens  3,  9,  and  10  possessed  some  slight  imperfections  in  bonding 
of  the  various  mosaic- type  Mocks  into  the  arch  configuration  and  accordingly 
exhibited  lower  strength  characteristics  than  did  better  constructed  speci¬ 
mens.  Specimens  1,  2,  4,  5,  6,  7,  and  11  collapsed  into  a  nonsymmetric 

configuration  upon  buckling.  For  all  of  these  specimens,  collapse  occurred 

2 

in  the  neighborhood  of  135  Ib/in  or  slightly  more,  which  is  in  good  agree¬ 
ment  with  theoretical  predictions  based  upon  an  asymmetric  buckling  mode. 

A  somewhat  greater  collapse  pressure  was  found  for  Specimen  8.  This  a*ch 
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2 

deflected  symmetrically  up  to  approximately  130  Ib/in  but  finally  buckled 
into  an  unsymmetric  configuration.  Figs.  13(c)  and  (d)  illustrate  this 
arch  as  progressively  greater  pressures  were  applied. 

The  experimentally  obtained  arch  displacements  were  usually  somewhat' 
greater  than  predicted  by  theory.  This  is  attributed  to  (i)  too  great 
a  resin  content  in  the  mosaic  blocks,  (ii)  nonhomogeneities  in  the 
blocks,  and  (iii)  a  fully  clamped  condition  was  not  achieved  at  the 
supports. 

An  overall  characteristic  of  the  arches  was  a  slight  brittleness. 

This  is  not  necessarily  an  inherent  property  of  such  a  composite  material 
but  is  believed  to  be  due  to  impurities  or  possible  dislocations  in  the 
specimens. 
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CHAPTER  V 
CONCLUSIONS 

The  analytical  results  obtained  from  the  classical “type  theory  lead 
to  higher  buckling  loads  if  transverse  shear  effects  are  included. 
Inclusion  of  transverse  normal  stress  effects  in  addition  to  transverse 
shear  still  further  decreases  the  predicted  buckling  loads.  The  predic¬ 
tions  of  these  three  types  of  theories  are  in  closer  agreement  as  the  thick¬ 
ness  of  the  arch  decreases. 

The  transverse  normal  stress  effect  on  a  thin  or  even  a  moderately 
thick  arch  is  small  and  can  usually  be  neglected.  However  the  effect  of 
transverse  shear  is  quite  significant  even  in  the  case  of  a  thin  arch. 

For  the  theory  incorporating  shear  deformation  effects  the  modification  in 
strain-displacement  relations  achieved  by  not  completely  neglecting  terms 
of  the  form  z/R  leads  to  predictions  in  better  agreement  with  experiments 
than  if  these  terms  are  completely  neglected. 

It  was  found  analytically  that  the  dimensionless  parameter  A  can  be 
used  to  characterize  regions  of  symmetric  and  asymmetric  buckling.  The 
value  of  10.8  for  this  parameter  separates  these  types  of  buckling  for  a 
hinged  end  arch  and  a  value  of  25  for  a  clamped  end  arch. 

Predictions  of  the  general  theory  developed  in  this  study,  which 
includes  transverse  shear  as  well  as  transverse  normal  stress  effects  are 
in  satisfactory  agreement  with  experimental  evidence  and  permit  reasonably 
accurate  predictions  of  load-deflection  relations  of  a  uniformly  loaded 
arch  fabricated  of  radially  reinforced  composite  material. 
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